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Abstract: 

We determine the strong coupling as from a next-to-leading order analysis of processes 
used for the NNPDF2.1 parton determination, which includes data from neutral and 
charged current deep-inelastic scattering, Drell-Yan and inclusive jet production. We find 
Us {Mz) = 0.1191 zbO.OOOG'^^P, where the uncertainty includes all statistical and systematic 
experimental uncertainties, but not purely theoretical uncertainties, which are expected 
to be rather larger. We study the dependence of the results on the dataset, by providing 
further determinations based respectively on deep-inelastic data only, and on HERA data 
only. The deep-inelastic fit gives the consistent result [Mz) = 0.1177 ± O.OOOD^^p, but 
the result of the HERA-only fit is only marginally consistent. We provide evidence that 
individual data subsets can have runaway directions due to poorly determined PDFs, thus 
suggesting that a global dataset is necessary for a reliable determination. 
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A precise knowledge of the value of the strong coupling constant Ug {Mz) [I] is nec- 
essary for accurate collider phenomenology, such as for instance Higgs searches at the 
Tevatron and the LHC In particular, in the gluon fusion channel the value of strong 
coupling is one of the dominant sources of uncertainty [3]. The current PDG [3] value 

a, (Mz) = 0.1184 ±0.0007 (1) 

is taken from Ref. [1], where it is obtained by combining several determinations, including 
some from processes (such as the r decay rate and the total e^e~ — >■ hadrons cross section) 
which do not require knowledge of nucleon structure, but others (such as deep-inelastic 
scattering) which do. The value of the uncertainty Eq. ([T|) may seem overly optimistic in 
view of the spread of values of available determinations and the significant dependence on 
the perturbative order of some of them: as a consequence, the use of a somewhat more 
conservative estimate of the uncertainty, such as Aa^ = 0.0012 at 68% confidence level 
has been recommended [2l[3] for LHC phenomenology. 

The determination of Ug {Mz) from the same wide of set of data which is used to 
determine PDFs is appealing because it simultaneously exploits the dependence on the 
coupling of scaling violations as well as that on individual hard matrix elements of the 
various processes under consideration. It is thus potentially quite accurate. On the other 
hand, in such a determination the value of is necessarily correlated to the best-fit 
form of the PDFs, and thus subject to potential sources of bias, such as for example an 
insufficiently flexible PDF parametrization. 

An example of the possible pitfalls of a simultaneous determination of PDFs and Ug 
is highlighted by the analysis of Ref. [7], in which the extraction of as from BCDMS and 
NMC deep-inelastic scattering data was performed using a methodology (scaling violations 
of truncated moments) which avoids completely the use of parton distributions. The result 
found, Us {Mz) = 0.124i!^Q QQ|, had rather different central value and uncertainties than 
those obtained by direct analysis of the same BCDMS {as {Mz) = 0.113 ± 0.005 [8]) and 
NMC {as {Mz) = 0.117l:o;oJg [9]) data by the respective collaborations. This suggests 
that the latter results, obtained using a PDF parametrization, were biased by it. 

Here we wish to provide a determination of exploiting the NNPDF methodology 
for determining PDFs |10ftl3| [T7] , which strives to avoid parametrization bias through the 
use of a Monte Carlo approach combined with neural networks as underlying unbiased 
interpolating functions. Specifically, we use the latest NNPDF set, NNPDF2.1 [18], which 
is based on a NLO global fit to all relevant hard scattering data, with heavy quark mass 
effects included through the so-called FONLL method [19]. NNPDF parton sets have been 
provided for a variety of values of as ■ Here we will use these sets to study the quality of the 
global agreement between theory and the data used in the PDF determination as the value 
of as is varied. A similar approach was used in Ref. [13] to provide a determination of the 
CKM matrix element |V^s| which turned out to be more accurate than any determination 
obtained from a single experiment. 

The use of NNPDF2.1 parton distributions has not only the advantage that parametriza- 
tion bias is reduced to a minimum, but also that the same methodology can be used to 
analyze different datasets, without having to retune the fitting procedure (such as, for 
instance, the form of parton parametrization) according to the size of the dataset. This 
enables a direct comparison of values of as obtained from different subsets of data which 
enter the global fit, and also an analysis of the correlation between individual datasets. 
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Figure 1: Left: The as a function of tlic number of replicas A^rop for NNPDF2.0; the horizontal 
line shows the value for A'rcp = 500. Right: The uncertainty a^2 Eq. ([2]) as a function of A^rop , 

— 1/2 

averaged over all values of ag . A fit of the form AN^cp is also shown. 



individual PDFs, and the value of Og- As a consequence, we will be able to address the 
issue of whether deep-inelastic scattering data systematically prefer lower values of as 
than hadron or e^e~ collider data. 

The main difficulty in determining as in the NNPDF approach is that, since PDFs 
are delivered in the form of a Monte Carlo sample, the quality of the fit (i.e. the of 
the comparison between data and theory) is a random variable, which only tends to a 
constant value in the limit in which the size of the Monte Carlo sample tends to infinity. 
The typical fluctuation of the for ^ single Monte Carlo replica is of the order of the 
square root of the number of data points A'dati while the fluctuations of the average over a 

1/2 

sample of A'rep replicas decrease as l/N^cp . So in order to be sensitive to variations of the 
total by a few units, as required for determination of a physical parameter, one needs 
for each value of as a number of replicas of the same order of magnitude as the number 
of independent data points. The total number of replicas required is thus rather large, 
which makes for a rather computationally intensive task. The value of the for & typical 
NNPDF fit (with A^dat = 3338) is shown as a function of TVrep in Fig. [H 

The uncertainty on the value of the due to the finite size of the replica sample may 
be computed using the so-called bootstrap method. Namely, the sample of N,-ep replicas is 
divided into A^part disjoint partitions with A'^rcp = -^rcpZ-^part replicas each. The variance 
of the x^ for tli6 full A'rep replica sample is then found from the variance of the A'part 
values x^ of each replica subsample according to 



2 _ 1 
x^) =W- 



part 




iVpart 



(2) 



The value of cr^2 , averaged over all the (eleven) values of Og to be considered, is displayed 
in Fig. [TJ A fit of the form ^A'rep , also shown in Fig. [H shows that the expected decrease 

1/2 

of the fluctuations with l/N^ep is borne out by the data. 

The determination of ag is performed by simply using a wide enough PDF replica set for 
Na^ fixed values of ag to compute for each value of ag the values of x^ and its uncertainty 
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Figure 2: The as a function of as(A/z) for tlie NNPDF2.1 global fit. The statistical uncertain- 
ties in the for each value of as have been determined from Eq. ([2]). The solid line is the result 
of a parabolic fit. 



(7^2 Eq. ([2]), and then fitting a parabola to the viewed as a function of Ug- The quality 
of the parabolic fit is then determined by evaluating the corresponding XparZ-^dofi with 
-^dof = — 3: a reasonable value of XparZ-^dof may be used to confirm that the parabolic 
approximation to X^(q^s) is adequate in the range of under investigation. The minimum 
of the parabola then provides the best-fit value of ag while the Ax^ = 1 range gives the 
uncertainty on it at a 68% confidence level. The further uncertainty due to the finite size 
of the replica sample is determined by error propagation of (t^2 Eq. ([2]) on the position of 
the minimum of the parabola. 

We now turn to results. For each value of we use a sample of at least A'rcp = 500 
replicas, with bigger samples of A^rep = 1000 replicas used for more sparse equally-spaced 
values in order to increase accuracy. For the fit to HERA data only, the range of values 
considered has been enlarged in order to ensure that the location of the minimum is 
approximately at the center of the region of ag which is being explored, and also because 
in this case the sensitivity to Ug is weaker due to the much smaller size of the data sample. 
The values of ag and numbers of replicas used in each case are summarized in Table [TJ 

The parabolic profile of \^ as a function of as{Mz) is shown in Fig.[2]for the NNPDF2.1 
global fit. Analogous results for fits to DIS data only and to HERA data only are shown 
in Fig. [31 The corresponding values and uncertainties of ag (Mz) are collected Table [21 In 
each case, we denote with "exp" the uncertainty from the Ax^ = 1 range and with "proc" 
the propagated "procedural" uncertainty, due the finite size of the replica sample. The 
quality of the parabolic fit is also shown in each case. 

The procedural uncertainties in Tab. [2l are all very small. This implies that the best-fit 
value of Og is already approximately independent of the size of the replica sample, and 
a further increase in the number of replicas is not necessary. However, it is useful to 
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Figure 3: Same as Fig. [2] but for a fit to DIS data only (left) and to HERA data only (right). 



2.1 global 





iVrop 


0.114 


500 


0.115 


500 


0.116 


1000 


0.117 


500 


0.118 


500 


0.119 


1000 


0.120 


500 


0.121 


500 


0.122 


1000 


0.123 


500 


0.124 


500 



2.1 DIS-only 



as (Mz) 


Afrop 


0.114 


500 


0.115 


500 


0.116 


1000 


0.117 


500 


0.118 


500 


0.119 


1000 


0.120 


500 


0.121 


500 


0.122 


1000 


0.123 


500 


0.124 


500 



2.1 HERA-only 



s (Mz) 




0.100 


1000 


0.102 


500 


0.104 


500 


0.106 


1000 


0.108 


500 


0.110 


500 


0.112 


1000 


0.114 


500 


0.115 


500 


0.116 


1000 


0.117 


500 


0.118 


500 


0.119 


1000 


0.120 


500 


0.121 


500 


0.122 


1000 


0.123 


500 


0.124 


500 



Table 1: The values of as (Mz) and the number of replicas iVrcp used in each case for various 
determinations of as {Mz). 



check this independence explicitly by varying the number of replicas, in order to make 
sure that the finite-size uncertainty has been determined correctly. To this purpose, we 
have repeated the three NNPDF2.1 determination of with A'rep = 500 for all values of 
as ■ Results are also collected in Table [2] and indeed show excellent stability: the change in 
value of Os is always smaller than the procedural uncertainty as the number of replicas is 
decreased. The values of the parabolic fit should follow a distribution with A'dof = 8 
degrees of freedom for the global and DIS fits, and A'dof = 12 for the HERA only fit. The 
standard deviation of Xpar/-^dof is thus expected to be of order 0.5, as indeed observed. 
We have finally checked that excluding the points at the edge of the fit, and adding extra 
parameters to the fit, has no significant effect on the results, and in particular it does not 
improve the quality of the parabolic fit. 

Our results for are displayed graphically in Fig. |H Our best-fit value of as is in 
good agreement with the PDG value Eq. ([1]), and has a surprisingly small experimental 
uncertainty. The experimental uncertainty increases as the size of the dataset is reduced, 
as it ought to. We see no evidence that DIS data prefer a significantly lower value of a^: 
the difference between values of the global and DIS-only determinations is of order of one 
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Table 2: Values of as {Mz) and associated uncertainties. All uncertainties shown are 68% confi- 
dence levels, with the experimental uncertainty obtained by requiring Ax^ = 1 about the minimum, 
and the procedural uncertainty from propagation of a^i Eq. ([5]) due to finite size of the replica 
sample. The quality of the parabolic fit as measured by Xpar/-^dof is also shown in each case. 
For the global, DIS-only and HERA-only fits (first three rows), the maximum number of replicas, 
given in Tab. [l] has been used. The three reduced replica fits (subsequent three rows) only differ 
from these because of the use of A^rop = 500 for all as values. The NNPDF2.0 fits of the last two 
rows also have A'j-op — 500 always. 

o", and thus entirely compatible with statistical fluctuations. Interestingly, the value found 
using HERA data only is much smaller, even though, because of the considerable (almost 
sixfold) increase in statistical uncertainty it is still less than three a from the global fit. 
The fact that HERA data prefer a lower value of ag may be related to the deviations 
between HERA data and the predicted NLO scaling violations which was observed in 
Refs. [5l|6] for the smallest x and HERA data: these may be affected by small x 
resummation or saturation effects. As shown there, scaling violations in this region are 
weaker than predicted from the behaviour observed in other kinematic regions, and thus 
would tend to bias the value of ag downwards. A dedicated analysis would be required to 
prove conclusively that this is the case. 

We make no attempt here to estimate theoretical uncertainties in our fit. Uncertainties 
due to inefficiencies of the global PDF fit (such as, for example, any residual bias related to 
parton distributions) should show up in the behaviour of the as a function of a^, either 
as point-to-point fluctuations or as a systematic deviation from the underlying unbiased 
quadratic behaviour (if they are correlated to the value of as). The good quality of the 
parabolic fit suggests that these uncertainties are small, and thus that our uncertainty 
is an accurate assessment of the total uncertainty due to the statistical and systematic 
uncertainties in the experimental data. On top of these, however, there will be genuine 
theoretical uncertainties related to the theory used in the computation of the various pro- 
cesses under investigation. Of these, the main ones are likely to be related to NNLO and 
higher QCD corrections (and possibly resummation of higher order QCD corrections in 
some kinematic regions), and to the treatment of heavy quark mass effects. These uncer- 
tainties are presumably of similar size here as in other determinations of as based on the 
same QCD processes; whereas they were studied systematically in older determinations 
(such as Ref. [I^), they have not been assessed for any of the more recent determinations. 
In cases in which both NLO and NNLO determinations are available, such as Ref. |15yi6j. 
a sizable downward shift of the best-fit value, of order of several percentage points, has 
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Figure 4: Left: Graphical representation of the values of as (Mz) of Tabic [2] (reduced replica fits 
not included). Uncertainties have been added in quadrature). Right: comparison to other recent 
determinations of as (Mz) from NLO PDF analysis. The PDG value of Ref. [1] is also shown. 




Figure 5: Same as Fig. [2] but for the NNPDF2.0 global fit (left) and NNPDF2.0 DIS only fit 
(right). 



been observed when going from NLO to NNLO. It will be interesting to see whether such 
an effect is also present in the NNPDF approach. 

As a very crude estimate of the order of magnitude of effects related to heavy quark 
masses, we have repeated the fit to the global dataset and that to DIS data using the 
NNPDF2.0 [T7] PDF set, which is based on a zero-mass variable flavour number scheme, 
in which all heavy quark masses are neglected. The results are also given in Table [2] and 
shown in Fig. [U while the corresponding parabolic fits are displayed in Fig [5l For these 
fits, Nj-cp = 500 replicas are used for all values of a^. Neglecting heavy quark mass effects 
induces a significant downward shift in a^. If one were to conservatively estimate the 
uncertainty due to heavy quark mass effects as the difference between the NNPDF2.1 and 
NNPDF2.0 results one would get, for the global fit, Aa^"^ fa 0.002. In the fit to DIS data 
only the shift is larger since the hadronic data are unaffected by the treatment of heavy 
quark mass effects. It is likely that the order of magnitude of uncertainties related to 
higher order corrections is comparable, so that the theoretical uncertainty is very likely to 
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Figure 6: The profiles for tire individual experiments in the NNPDF2.1 global fit together with 
the results of the corresponding parabolic fits to Qg. The uncertainties due to the finite size of the 
replica sample are shown on each value. 



be the dominant one. 

Our result using the standard (global) NNPDF2.1 set is compared in Fig. [5] to other 
recent NLO determinations of Os which rely respectively on the MSTW [H], CTEQ [21] 
and ABKM [16] PDF sets (ABKM and MSTW also provide determinations at NNLO). The 
MSTW and ABKM groups perform a simultaneous fit of PDFs and the strong coupling, 
thus obtaining a correlated Hessian matrix which mixes the PDF parameters with , while 
CTEQ simply studies the dependence of the fit quality on as is done here (see Sect. 3 of 
Ref . [20] ) . The equivalence of the q ^ uncertainty obtained from either method is explicitly 
shown in Ref. [21] . The dataset on which the CTEQ and MSTW determinations are based 
is very similar to our own, while ABKM use a smaller dataset, which in particular does 
not include collider jet and vector boson production data. All these determinations are in 
agreement with each other within uncertainties. The rather larger statistical uncertainties 
found by CTEQ and MSTW can be understood as a consequence of the fact that these 
groups [151121] use a tolerance [22] criterion to obtain 68% confidence levels, based on a 
substantial rescaling of the uncertainty ranges in parameter space. This is not necessary 
in our approach because, once the is treated as a random variable, its fluctuations can 
be studied (and in particular kept under control) by a suitable choice of the size of the 
Monte Carlo sample, as discussed above. 

This conclusion may be cross-checked using a variant of the method suggested in 
Ref. [23j , namely, by checking whether the distribution of results obtained from individual 
datasets follows a gaussian distribution and determining the width of this distribution. 
To this purpose, we have performed a parabolic fit to the profile for each experiment 
entering in the global NNPDF2.1 determination. Results are displayed in Fig. [6] together 
with the uncertainties due to the finite size of the replica sample, determined as above. 
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Distribution of pulls with global fit a3~| 




Figure 7: Distribution of pulls, Eq. ([3]), for the value of preferred by the individual experiments 
included in the global fit. These pulls have been summarized in Table 15]). 



These uncertainties are clearly much larger than the point-to-point fluctuations of the 
individual values for each experiment, due to the fact that the latter, being determined 
from a global fit, are strongly correlated with each other. The value of and its statistical 
uncertainty for each experiment are determined by performing a parabolic fit to each 
profile: results are collected in Table El for all experiments for which there is a minimum 
in the fitted range. This is not the case for the NMCratio, SLAC, CHORUS, H1F2C and 
FLH108 data. 

The distribution of results can be studied defining the pull 

^ ai(Mz)-ar (Mz) ^ 

/ i,2 , tot, 2 
Y (Jas + 

where a* (Mz) is the best fit value for the i~th experiment and the associated statistical 
uncertainty, obtained from the Ax^ = 1 rule. The pulls are summarized in Table [3] and 
displayed graphically in Fig. [71 A gaussian fit to the distribution of pulls is performed, 
and also displayed in Fig. [71 The gaussian fit is in good agreement with the histogram 
data with mean (P) = 0.04 and standard deviation ap = 1.3. The standard deviation 
would be further reduced somewhat if finite-size uncertainties were included; this however 
would require a lengthy correlation analysis. We conclude that the value of the tolerance 
required to get a perfectly gaussian distribution of pulls is smaller than 1.3 — a value 
which is clearly compatible with a statistical fluctuation. 

Finally, we exploit the fact that in our approach a single procedure can be used to 
obtain PDFs from datasets of different size, in order to study the issue of whether (and 
why) different values of may be preferred by different datasets. In Fig. [8lwe compare the 

profiles for the global NNPDF2.1 fit, already shown in Fig. [6l to the same quantities 
determined for the fit to DIS data only. The behaviour of the fit quality for BCDMS 
data is particularly interesting: these data have been repeatedly found [8l[T5l|23] to prefer 



9 



Experiment 
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P. 


NMCp 


0.1192 ±0.0018 


-0.05 
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0.1204 ±0.0015 


-0.78 


HERA-I 


0.1223 ±0.0018 


-1.65 
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0.1170 ±0.0027 
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0.1252 ±0.0068 
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0.1168 ±0.0100 


0.22 
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0.1135 ±0.0029 


1.87 


CDFWASY 


0.1181 ±0.006 


0.16 


CDFZRAP 


0.1150 ±0.0034 


1.18 


DOZRAP 


0.1227 ±0.0067 


-0.53 


CDFR2KT 


0.1228 ±0.0021 


-1.67 


D0R2CON 


0.1141 ±0.0031 


1.57 



Table 3: The pulls Pi Eq. ([3]) for individual experiment included in the NNPDF2.1 global fit case, 
computed for each experiment which has a minimum in the range considered. 

a relatively low value of as (in particular, lower than Eq. ([ID). It turns out that the 
profile for these data is rather different according to whether one is looking at PDFs 
determined using DIS data only (green, lower curve in Fig. [8]) or a global fit (red, higher 
curve in in Fig. [8]). Indeed, in the DIS-only case these data indeed seem to prefer a lower 
value of a<j, but this is no longer the case in the global fit. This suggests that as is 
lowered, the quality of the fit to BCDMS data can be improved by changing the PDFs in 
a way which is allowed by DIS data, but which is not compatible with other data in the 
global fit. In other words, there is a runaway direction in the space of PDF parameters 
along which the for BCDMS data decreases as is lowered, which is forbidden in the 
global fit because then the for some other dataset would increase. 

The situation can be further elucidated by studying the correlation [18] between par- 
ton distributions and the value of the for individual experiments. The presence of a 
nonvanishing correlation means that, at the best fit, the x^ for that experiment is not 
stationary, i.e. it can be lowered or raised by changing the given PDF. Correlations of 
opposite sign for different experiments then mean that these experiments are pulling the 
PDF in opposite directions. The correlation coefficients are shown for the gluon PDF as 
a function of in Fig. [9] for a pair of values of x, for the global NNPDF2.1 PDF set (all 
computed from a set of A'l-cp = 500 replicas). It is apparent that while for larger values 
of Os ~ 0.120 correlations are small and with the same sign, as ag is lowered correlations 
become larger, with opposite sign for jet and DIS (HERA and BCDMS) experiments (with 
the Drell-Yan experiment E866 showing no significant correlation). This means that in- 
deed, as suggested above, for low as DIS and jet experiments pull the gluon in opposite 
directions, while they become more consistent for larger a^: hence, a determination of 
including DIS data only can easily be biased. The fact that BCDMS data prefer a lower 
value of Os in a DIS-only fit, but not if the gluon is constrained by jet data was also found 
recently in Ref. [2l], in the context of the MSTW08 parton determination. However, in 
that case the BCDMS data were also found to significantly bias downwards the value of 
as of the DIS fit, perhaps due to the fact that the MSTW gluon parametrization, though 
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Figure 8: Comparison of the profiles for the global NNPDF2.1 fit (same as in Fig. [51 red, solid 
curves) to those determined for the DIS-only NNPDF2.1 fit (green, dashed curves). 



more flexible than that of other groups, is still less flexible than that of present analysis. 

The fact that runaway directions for the may appear in the joint a^-gluon space 
can be understood by noting that in DIS the gluon is determined by scaling violations, 
hence a smaller value of Og can be partially compensated by a larger gluon and conversely. 
However, the jet cross section pins down the size of the gluon (at the rather larger scale of 
the jet data) thereby quenching this potential instability. Hence, we conclude that even 
though in our fit the DIS-only value of is not significantly smaller than that for the 
global fit (possibly due to the great flexibility of the functional form of our PDFs), a fit 
to DIS data, and specifically to BCDMS data, has a potential instability in the direction 
of lower values of Og which is only kept under control by the inclusion of jet data. 

We conclude that a reliable as determination, with surprisingly small statistical uncer- 
tainty, can be obtained by a combined analysis of a wide set of data which simultaneously 
depend on the value of the strong coupling and the parton distributions. Theoretical un- 
certainties are likely to be dominant and significant. They could be kept under control at 
least in part by inclusion of higher order corrections, in particular NNLO as well as all- 
order resummation, which might be especially relevant for small- 2: HERA data [5l[6tl25j. 
Once resummation corrections are properly included, it might be convenient to use for 
DIS only HERA data, which are free of ambiguity related to nuclear corrections or power- 
suppressed corrections. It will be interesting to repeat the analysis presented here once 
PDF sets which include these higher order effects become available. 
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Figure 9: Correlation between tlie and the input gluon as a function of ag (Mz) for x = 0.05 
(left) and x = 0.1 (right) for the NNPDF2.1 global PDF set. 
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